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Abstract 

We continue to investigate the minimal 5*0(10) x U{1)h model which we pro- 
posed recently. Renormalization group analysis of the model results in natural 
predictions of quark-lepton masses and Kobayashi-Maskawa matrix along with 
neutrino mixings adequate for solar neutrino oscillation. 
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1 Introduction 



Grand unification has been a fascinating idea for two decades. Among the possible gauge 
groups, 5*0(10) is the smallest candidate that can incorporate the observed fermions of one 
generation into an irreducible multiplet: it attains matter unification for a single generation 
of quarks and leptons. However, SO{10) grand unification by itself provides no natural place 
for triplicity of generations, to say nothing of the hierarchical structure of mass matrices. For 
instance, it gives no explanation on the fact that the mass of top quark is more than 10^ times 
that of electron. 

In a previous paper we constructed an 50(10) x U{1)h model as an attempt to place 
generation structure on a plausible position in grand unification. We introduced a minimal 
Higgs content to break the gauge symmetry without any additional scalars. The horizontal 
Peccei-Quinn symmetry U{1)h distinguishes three generations of fermions to impose Georgi- 
Jarlskog(GJ) relations at the unification scale M[/. It also makes Yukawa couplings of 
order one be able to realize hierarchical mass matrices with the aid of remnant effects of 
certain irrelevant terms suppressed by the Planck scale Mp as the cut-off scale in the theory 
(For more details, see ref. [0]). 

The model has three factors which determine the values of mass matrices in combination 
with order-one Yukawa couplings at the unification scale. The first one is the hierarchy 
factor e = Mu/Mp brought about by the remnant effects mentioned above. The second one 
comes from the mixing of Higgs bosons which is represented by parameters a, /3, 7, 5 defined 
later. The third is the running of Yukawa couplings according to renormalization group (RG) 
equations. 

In this paper, we proceed to analyze mainly the third factor to compare predictions of our 
model with experimental estimates. Numerical analysis of RG equations results in natural 
predictions of quark-lepton masses and Kobayashi-Maskawa(KM) matrix along with neutrino 
mixings adequate for solar neutrino oscillation 

The paper is organized as follows: We first recapitulate the setup of our model in section 
2 to derive RG equations in section 3. We restrict ourselves to dealing with one- loop RG 
equations and ignore threshold corrections throughout the paper. Section 4 makes exposition 
of numerical results obtained by RG analysis. Section 5 concludes the paper. Some definitions 
are given in Appendix A and analytical consideration on quark mass matrices is made in 
Appendix B. 
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2 Yukawa Interactions 

In this section, we present effective Yukawa interactions under the chain of symmetry breaking 
considered in ref. . There are three scales of symmetry breaking postulated (with the help of 
fine-tuning) in the model: the unification scale Mjj, the intermediate scale M/, and the weak 
scale My/. 

We introduced Higgs fields $(210,-8), A(126, — 10), and iJ(10, — 2) transforming under 
50(10) X U{1)h- The field $(210, -8) develops a VEV of order Mu, which breaks 50(10) x 
U{1)h into G224 X Zg, where G224 denotes the Pati-Salam group SU{2)l x SU{2)r x SU{4)c. 
In particular, D-parity oddness of $(210, —8) causes violation of the left-right symmetry. 

We assumed that there remain iJ(2,2,l), A(2,2, 15), and A(l,3,10) on the Pati-Salam 
stage between the scales Mjj and M/, the other Higgs fields with masses of order Mu de- 
coupling from the system. This is a modification of the minimal fine-tuning, which claims 
non- decoupling of only if (2, 2, 1) and A(l, 3, 10). The additional field A(2, 2, 15) develops an 
induced VEV of a considerable size to contribute fermion mass matrices [^. 

The effective Yukawa interactions on this stage are given by 

+Y3l/jl(pl5^R + n^I^lsV'R + Y5^l(j)io^R + h.c, (2.1) 

where F's are Yukawa couplings and summation over suppressed generation indices should be 
understood. Fermions ipL and ipR denote (2,1,4) and (1,2,4) representations, respectively. 
Scalars 0i, 0i5, and 0io correspond to H{2, 2, 1), A(2, 2, 15), and A(l, 3, 10) in this turn (See 
Appendix A). 

Yukawa couplings at the unification scale Mjj provide a boundary condition 

/ -e^yn £^13 \ / 6^yi2 \ 

Yi=\ 0, ^2= e^yi2 -5^23 , 

V -eyi3 / \ ey23 / 



/ Zi3 \ ^ / ezn 2:13 \ 

Y3=\ Z22 \, n = 0, ^5 = ^ ^22 , (2.2) 

V ^13 ; V Z,3 eZ33 ) 

where Y's have been regarded as matrices in generation space. The small factor e = Mjj/Mp 
stems from the remnant effects and y's and z's are input parameters of order one, which we 
expect to be determined by more fundamental theory presumably including gravitation. GJ 
relations can be obtained from the above texture as shown in ref. [0 (See also Appendix B). 
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Below the intermediate scale Mj, the model effectively coincides with the standard model, 
and the effective Yukawa interactions are given by 

£2 = (Q 0) dR + (Q 0) + Y} (L(t>) cr + ur + h.c, (2.3) 

where Q, and L denote the standard Higgs, quark, and lepton doublets, respectively. 
Although this Lagrangian contains Dirac mass terms for neutrinos Y}{L(I))vr, they can be 
approximately neglected below Mj due to order-M/ Majorana masses of right-handed 
neutrinos vr, which originates from the term Y^ip'^i^iQipR in eq.( p.l|) . 

The standard Higgs doublet is a linear combination of four doublets contained in 
i/(2,2,l) and A(2,2,15): 



= ai/i + +7A1 +5A_i (2.4) 

2222 

with a normalization condition 

a' + /32 + 7' + 52 = 1, (2.5) 

where subscripts ±| indicate hypercharges of the fields. Hence the following relations hold at 
the intermediate scale Mf. 

y„ = aY,+pY, + ^^Y, + ^5Y,, 

Ya = aY2 + PY, + ^^Y, + ^SYs, 

= aY, + /3Y2 + 0=^Ys + ^SY,, 

Y, = aY2 + PY^ + ^^Y, + ^SY,, (2.6) 

where the coefficients and result from the normalization of 0i5 defined in the Appendix 
A. The parameters a, (3, 7, 5 are written in terms of the VEVs defined in ref.[|I| as follows: 

a = ^, /? = ^, 7 = 2^3^, 5 = 2v^^, (2.7) 

where vw denote the VEV of the standard Higgs doublet. Note that the VEVs w* and w* are 
of order e relative to vw under the assumption in ref.|jl| that they are induced in the Higgs 
015 with mass of order Mj. Hence 7 and 6 take values of order e. This makes it natural to 
define order-one quantity 7' and 6' in terms of 

7 = 7'e, 6 = 6'e. (2.8) 
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It should be emphasized that only two parameters e and (3 are small to realize the hierar- 
chical structure of mass matrices (See Appendix B). The other parameters y's, z's, 7', and 5' 
are of order one, which merely affect detailed numerical values of masses and mixings without 
altering their orders of magnitude. 



3 Renormalization Group Equations 

In this section, we show one-loop RG equations for gauge and Yukawa couplings ((^'s and 
y's) in the effective theories described in the previous section. 

We first give RG equations for the 6*224 theory on the stage between the unification and 
intermediate scales: 

"^^-^W, ''<=(2.f,4). . = 2..2h,4c, (3.1) 

where cUj = t = In/i, and /i denotes a renormalization point in the MS scheme; and 



dYi 



dY 

IGtt^^ = Y2PL + PiY2 + Y2Pi + I3,2 + Y2[3g, 
dY 

dY 

dY 

levr^^ = Y,(3R + (3iY, + Y,(3io + Y,(3'g, (3.2) 

where (3li I3r and 15, 10 correspond to contributions from wave-function renormalization of 
V'L, "^i?, and 01, 15, 10, respectively; /?„'s and [3'g correspond to contributions from vertex 
renormalization and gauge couplings: 

(5l = yIy^ + Y^Y2 + ^{Y^Y, + yIy^), 

(3r = Y^Yi + Y^Y2 + ^iYiY3 + Y^Y^ + Y,%), 

A = 4tT{Y,^Y, + Y^Y2), 

I3i, = ti{Y^Y, + YlY,), 

/3io = iT{Y^,Y,), 

(3,, = -2Y,Y^Y2-2Y2Y^Y,~^Y3Y^Y^-^Y,Y^Y,, 
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15, 

y 



-2Y2Y^Y, - 2Fiy/F2 - ^YsYIy, - ^Y.Y^Y,, 



\y^yIy^ + ^41^4^3 - 2Y^YIY2 - 2Y^YIY^, 
]^^YlY^ + ]y^yIy^ - 2Y^YlY^ - 2Y^Y^Y^, 

9 2 2 



(3.3) 



Let us turn to the energy region below the intermediate scale. RG equations on this stage 
are those for the standard model: 



duji 



27r 



41 

To' 



19 



-7 



3c; 



(3.4) 



IGtt 
IGvr 
167r' 



2dYd 



dt 
,dYe 
dt 



Yd 



(riK - Y^Ya) + tr(3FiF. + 3^]^, + r;!;; 



17 
20 



9i 
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+ 8^/1 



^-{YlYd - Y^Y^) + tr(3y„ty„ + W^Y^ + F^Fe) - (^-g^ + + 8^73^ 

^-yX + tr(3y;ty„ + srjr, + y}y,) - ^(^^^^ + gl 



For a recent analysis of them, see Ref . [ 



4 Numerical Results 



(3.5) 



We proceed to consider numerical solutions to the RG equations listed in the previous section. 
The running of the gauge couplings are independent of Yukawa couplings to the extent of one- 
loop analysis. Thus we can first obtain the unification scale Mu and the intermediate scale 
by means of RG equations for gauge couplings (3J) and (3^) with their experimental 



values below the weak scale 110] as a boundary condition: 



Mr, 



10^^-^GeV, Mr 



10 



11.2 



GeV. 



(4.1) 



The unified gauge coupling at Mjj comes out to be gu ~ 0.585. We think of this as a typical 
value of order one, since gauge coupling seems fundamental in view of its geometrical origin. 
Yukawa couplings at the unification scale are to be compared with this value as the standard 
one. (Conversely, one might also regard this value as an experimental evidence of coupling 
constants being of order one.) 
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Now that we have obtained running gauge couphngs, we turn to analyzing Yukawa sector. 
We assume CP conservation in the following analysis. In particular, KM phase is set to 
zero. Thus the case of large KM phase are excluded from the analysis in this paper, though 
small phase may be taken into account perturbatively and seems not to affect the results 
considerably. 

The procedure we pursue is as follows: To begin with, we choose appropriate values for the 
input parameters e, /3, j/'s, z's, 7', and 5' partly with the help of trial and error (See Appendix 
B). Note that y's, 2;'s, 7', and 5' must be of order one. Then we make Yukawa couplings evolve 
from the unification scale Mjj down to the weak scale M\y ~ 174GeV by solving RG equations 
numerically. Finally the resultant Yukawa matrices at the weak scale are diagonalized to yield 
fermion masses and mixings. Neutrino masses and mixings are derived by means of sea-saw 
approximation Q from the values of Yukawa couplings at the intermediate scale M/, where the 
right-handed neutrinos are supposed to decouple. We compare the results with experimental 



estimates of running masses and mixings at the weak scale |TT |. 

Let us exhibit a numerical sample which provides a realistic pattern of mass matrices at 
the weak scale: input parameters in table lead to the results in tables |^ and ^. The masses 
of neutrinos are given by 
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K„ m. , m,J ~ (4.1 X 5.9 x 10-^ 3.0 x lO"') x ^, (4.2) 

Vi 

where vi denotes the VEV developed by the Higgs 0io. The values vw — 174GeV and 
vi ~ lO^^-^GeV predict 

m,^ ~ 10-^, (4.3) 

and a negligible value of tjIu^ relative to m;^^. This is consistent with the small-angle MSW 
solution to the solar neutrino problem 0, which is implemented by 

Ame^ = (m^^ - m^J 1/2 ^io-3eV, 

sin^e/. ^ 0.03 -0.06. (4.4) 

The above example shows that qualitative agreement is achieved between predictions of the 
model and experimental estimates. In particular, the hierarchical structure of mass matrices 
was shown to be indeed realized in terms of Yukawa couplings exclusively of order one. 

So far so good. However, we cannot help complaining about the prediction of bottom mass. 
It has tendency to come out larger [|12| in this model than the experimental estimate (provided 



the tau mass is fitted.) This defect originates from GJ relation mf, ~ at the unification 
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1^1 = 



/ 0.3 X £2 0.3 X £ 


\ -0.3 X £ 0.5 



( 3.25 X £2 \ 
3.25 X £2 0.16 X £ 
-0.16 X £ 



1 

250' 



/ 0.5 X £ -0.145 \ 

0.65 
V -0.145 0.5 X £ y 



/3 



53' 



7 



25 
"9"^ 



5' 



Table 1: Sample input parameters at = 10^^'^GeV 





Prediction 


Experiment 




Prediction 


Experiment 


mu 


2.5 X 10-3 


2.4 X lO--'^ 




3.3 X 10-3 


4.2 X 10-3 


rUc 


0.61 


0.61 


ITLs 


0.087 


0.085 




160 






3.7 


2.9 




Prediction 


Experiment 




Prediction 


Experiment 


me 


5 X 10-* 


5 X 10-4 




~ 4 X 10-^9 






0.1 


0.1 




~ 6 X 10-^3 




rrir 


1.7 


1.7 




~ 3 X 10"^° 





Table 2: Running quark and lepton masses (GeV) at = 174GeV. 



quark 



lepton 



-0.98 
-0.22 
-0.005 

1 

0.059 
0.014 



0.22 
-0.98 
-0.052 

-0.058 
1 

-0.063 




Table 3: Predictions of quark and lepton mixing matrices aX ji = 174GeV. 
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scale and is scarcely dependent on the parameters we can choose to get realistic predictions. 
Taking this issue seriously, we might need some modifications of the present model. Even in 
such circumstances, we hope that the qualitative features of the model survive to naturally 
achieve the mass-matrix hierarchy. 



5 Concluding Remarks 



In the preceding sections, we have investigated the running of Yukawa couplings below the 
unification scale in our model of grand unification The remnant effects of orders e and 
play crucial roles to make realistic predictions in the model. Without the remnants, the model 
would possess a parity symmetry for the second generation iIj2 —>■ —'4'2i which completely 
forbids mixings between the second and the other generations: that is, in the renormalizable 
setting, radiative corrections could not generate necessary operators corresponding to the 
remnants. 

Finally let us see the behavior of couplings above the unification scale. The Yukawa 
interactions are given by 

where higher-dimensional terms are ignored. Greek indices are 5*0(10) vector ones, 7^ yields 
32-dimensional representation of Clifford algebra { 7^ , 7^ } = 2 5^,^, and B denotes charge 
conjugation for S'O(IO) spinors: B = Yi/i-.odd 7/^- 126 Higgs A is self-dual antisymmetric 
tensor, which satisfy 

_ 'i _ 

where ^^^^2■■■^lw denotes the invariant antisymmetric tensor. 

The one-loop RG equation for the unified gauge coupling g is given by 

which shows that the theory is asymptotically non-free. Note that the absolute value of the 
beta function is so small that the coupling stays within perturbative range up to the Planck 
scale. The running of effective gauge couplings are exhibited in figure |T]. 

The Yukawa couplings above the unification scale obey approximate RG equations 

IGtT^ j ° = -— ^10^10^10 + -TT ^126^126^10 + ^10^126^126 

at IQ 8 8 

+Fiotr(rioFil,)- 24(7^^10, (5.4) 
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16vr^ — — — — ^126^126^126 + ^10^10^126 + T77 ^126^10^10 

at 4 Id Id 

+2ri26tr(ri26n26) - 24(7'rio (5.5) 

At the unification scale, they satisfy 

(^10)33 = (1^1)33, 2^126 = ^3. (5.6) 

The flow of Yukawa coupling (1^10)33 with respect to the unified gauge coupling g is shown 
in figure 0. The Plank scale corresponds to ~ 0.62 and the GUT scale to (7 ^ 0.585. Thus 
the Yukawa coupling 7/33 in the numerical sample in the previous section is found to be of 
order one even at the Planck scale. This is consistent with the perturbative treatment above 
and the general philosophy of effective field theory that coupling constants are of order one at 
the cut-off scale. 

In fact, figure ^ suggests that (1^0)33 is almost always of order one at the GUT scale 
whatever it is at the Plank scale. One can even consider the case where the Yukawa coupling 
blows up at the Plank scale, which might indicate the presence of some dynamical phenomenon 
out of perturbative picture. 
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0.56 0.58 0.6 0.62 0.64 0.66 0.68 0.7 

9 



Figure 2: Running Yukawa coupling above Mu\ 
g denotes the unified gauge coupling. 
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Appendix 



A Definitions 

On the 6*224 stage, the following representations of the fields under SU{2)ji x SU{2)l are 
employed: 

01,15 fo^ , 0io = 0roT"^; r^^^a^, 1,2,3 (A.l) 

where U and D represent quartets of SU{4)c and a"^ denote Pauli matrices. The definition 
of 01,15 is given by 

01,15 = e'^0M5e; ^ = ( -1 ) ■ ^^-^^ 
The SU{4:)c representations of 0i5 and 0io are given by 

15 15 ig 15 1 

a=l a=l ^ 6=1 ^ 

10 10 p- 

010 = E </'?o'5a; tr(5<,5'^) = C E'5a-5"=ol' S'^ = Sl, (A.3) 

a=l a=l ^ 

where T" correspond to 5't/(4)c generators in the defining representation. 

B Quark Mass Matrices 

It seems instructive to consider some analytical relations among quark masses and mixings 
which are expected to hold at the weak scale My/ in our model. We make a rough estimate 
that the running of each Yukawa coupling between the scales Mjj and My/ does not affect its 
order of magnitude, which can be checked numerically. Then the mass matrices at the weak 
scale can be written in terms of rescaled couplings 

1/33 ~ ^33 ~ «2/33, ^22 ~ 4 ~ ^^^2, etc. (B.l) 
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as follows: 



M,, 



ez. 



22 



e^y23 

\<^'l3-yi3) ^iy23 1/33 

e'rVu ^r'C4 -^r'yls 



where 



I c(rX43 - vis) ^rvi 



a y 



(B.2) 



(B.3) 



Eq.( ^.5| ) implies that the ratio between the Yukawa couplings of top and bottom quarks does 
not change considerably 

Z/33 ^ 1/33 (B.4) 

under dominance of the gauge couplings g2 and g^. 

We now derive approximate relations for the parameters ^ and C with the aid of smallness 
of the hierarchy factor e = Mu /Mp ~ 1/250 (See section 3). The masses of top and bottom 
quarks are expressed as 

^ = l/33 + 0(e^), ^-^-^ ' --^-^^ 

which implies 



1/33 + 0{e'C 



(B.5) 



C ^ — • 

nit 

If we put rrit ~ 160GeV, we have C, — 1/50, which yields e^^^ ~ 1/5. 
Let us proceed to the second generation. We obtain 



(B.6) 



_2fl/23^ 



Vw yh 



er\42 + e\-^y^ 
1/33 



3t-3\ 



with a KM matrix element 



Vr 



cb 



1/33 



.y^3^ 



yh' ^ 



(B.7) 



(B.8) 



Note that the ©(e^^-^) term in eq. (|B.7|) is retained in anticipation of smallness of (. By means 
of these relations, we get 



rrir 



(B.9) 



which is satisfied when ~ 1/10. The enhancement factor ^"^ for the mixing Vcb in eq.( p.8|) 
comes from contribution of H, whose coupling is characteristic of non-supersymmetric models. 
A choice yi^^/y^s ~ 1/4 yields Vcb - 1/20. 
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